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Abstract: In this paper, we propose a novel stacked laser dielectric acceleration structure. This 
structure is based on the inverse Cherenkov effect and represented by a parametric design 
formulation. Compared to existing dielectric laser accelerators relying on the inverse Smith- 
Purcell effect, the proposed structure provides an extended-duration synchronous acceleration 
field without requiring the pulse front tilting technique. This advantage significantly reduces 
the required pulse duration. In addition, the easy-to-integrate layered structure facilitates 
cascade acceleration, and simulations have shown that low-energy electron beams can be 
cascaded through high gradients over extended distances. These practical advantages 
demonstrate the potential of this new structure for future chip accelerators. 


1. Introduction 

The industrial applications of modern particle accelerators have been extended beyond 
scientific research [1,2]. However, radio-frequency particle accelerators are typically bulky and 
heavy with limited practical limitations, such as the threshold for metallic structure breakdown 
and restrictions in the power and wavelength of microwave sources [3-6]. In recent decades, 
researchers have attempted to construct miniature and even benchtop particle accelerators [7- 
15]. Acceleration mechanisms based on laser-dielectric interactions have recently attracted 
increased attention, with their promising potential to achieve the development of on-chip 
microscopic accelerators. The inverse Smith—Purcell-effect dielectric laser accelerator (ISP- 
DLA) [8-15] and the inverse Cherenkov-effect dielectric laser accelerator (ICS-DLA) [16-19] 
are two laser-dielectric acceleration concepts that have been investigated since the 1960s, 
shortly after the invention of the laser. 

Most dielectric laser accelerator (DLA) schemes are based on the inverse Smith—Purcell 
effect. In their typical configuration, a periodic dielectric structure is adopted to spatially 
modulate the incident laser field and synchronize the incident-free electrons. This is achieved 
through the spatial harmonics of the induced waves in the electron beam channel. These 
previous developments have several challenges [8-15], including the following: (1) microscale 
fabrication of accelerator structures, which requires advanced micro-processing and nano- 
processing techniques; (2) selection of dielectric materials, some of which have high laser 
damage thresholds but cannot be fabricated as microstructures and nanostructures; (3) reduced 
laser damage thresholds of dielectric materials after nano-processing; (4) the pulse front tilting 
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technique, which is required to extend the action time of the laser and dielectric to achieve 
extended acceleration, improve laser energy utilization efficiency, and prevent unnecessary 
high energy loads on the structure [20]. These challenges limit the acceleration gradient, 
acceleration distance, energy utilization efficiency, and flexibility of laser dielectric 
acceleration structures while increasing experimental difficulty. 

In this study, a stacked acceleration structure based on the inverse Cherenkov effect was 
developed to achieve cascade acceleration of nonrelativistic low-energy electrons through a 
multilayer dielectric. First, the proposed structure is identical to the existing inverse Cherenkov- 
effect accelerator and provides effectively solves the challenges encountered in the acceleration 
structure of the inverse Smith—Purcell effect, including (1) elimination of the requirement for 
laser frontier technology and (2) not requiring multiple lasers to achieve cascade acceleration, 
eliminating the difficulty of coupling multiple lasers. Second, this structure considers higher 
acceleration gradients and energy gains and has no theoretical upper limit on acceleration 
energy. In addition, because this stackable structure can be assembled from modular 
acceleration units, it provides additional configuration flexibility to satisfy specific 
requirements for different applications. This flexibility of the stackable design makes the 
current structure a prototype for future chip accelerators in universal application scenarios. 


2. Theory and structural design 

The proposed stacked acceleration structure is depicted in Fig. 1. A case study of two 
dielectric layers is presented in this paper for demonstration without loss of generality. The 
line-shaped collimated incident laser was perpendicular to Surface II and the beveled surface 
of a multilayer dielectric prism with a right-angled triangular cross-section. An evanescent 
wave is generated at Surface I with a velocity of v, = c/(n: sin æ), travelling forward along 
the surface in the x-direction, where æ is the angle between the hypotenuse and the bottom 
edge, c is the velocity of light in vacuum, and n is the refractive index of the prism. This 
wave has a longitudinal electric field component E,, transverse electric field component Ey, 
and out-of-plane magnetic field component Bz. 
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Fig. 1. (a) Schematic of stacked acceleration structure. The red arrow line represents the incident laser. The 
refractive indices of the two dielectric layers are n, and nj. (b) Partial enlargement of Fig. 1(a) showing refraction 
process occurring at interface. The angles between the bottom edge of the right-angled triangle and the two traversal 
directions perpendicular to two laser beams (equivalently, the electric field directions) are œ and az, respectively. 


The electrons move in the x-direction in the vacuum below Surface I. Electrons in phase 
with an accelerating electric field (equivalently, negative phase Ey, <0) can then be 
accelerated by applying the inverse Cherenkov effect. However, the electron velocity v, must 
be equal to that of the evanescent wave, that is, Ve = Vp, to achieve continuous acceleration. 
On Surface I, the first criterion [17-19] for achieving continuous acceleration through the 
inverse Cherenkov effect is 


sina =c/(ven) = 1/(Bn). (1) 
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In Equation (1), æ is the angle between the hypotenuse and the bottom edge, c and ve 
are the velocities of the light in vacuum and electrons, respectively, and B = v/c is the 
relativistic velocity of the electron. As the electron accelerates and v, increases, the electrons 
move faster than the wave, lose synchronization, and eventually, their acceleration end. 
Currently, the solution to this problem is achieved by varying angle a [18-19]. In this study, 
the electron velocity was matched to the wave velocity by selecting a different dielectric to vary 
the refractive index. First, the refracted light ray A is analyzed, as shown in Fig. 1. Each time 
the laser passes through the interface between two adjacent layers with different dielectrics, the 
light is refracted owing to a change in the refractive index. Some of the laser light is reflected, 
which reduces the laser energy; this reduction is undesirable for the energy utilization efficiency 
and accelerating electric field. Thus, it is necessary to ensure that the light passes through the 
interface only once so that the laser experiences only one additional reflection and refraction. 
A limit was proposed for the dielectric length of each layer to achieve this. Taking the second 
layer in Fig. 1 as an example, the second-layer thickness Lz is limited by Rays B and C. The 
geometric relationship is expressed by L, = L, X (tan a)”. Similarly, the recurrence equation 
for the thickness L; of each subsequent layer can be derived, Lj4, = L;/(cosa@)*,i => 2. By 
combining these two equations, the relationship between the thickness L; of each layer and 
the thickness L, of the first layer is expressed by L; = L4 x (tana)*/(cos a)?“ ,i > 2. 
Second, even when the laser passes through the interface only once, as shown in Fig. 1, the 
angle of the laser to the bottom surface (Surface I) changes from œ to æ, after refraction at 
the dielectric interface. Equation (1) must still be satisfied to maintain acceleration: sina, = 
1/(B,nz). When refraction occurs at the interface, the law of refraction is expressed by 
n,sin@, =n, sin, . In addition, the trigonometric relationship is 1 = (sina@,)? + 
(cos a3)? = [1/(B2n2)]? + (sin 02)? = [1/(B2n2)]? + (n; sin 01/n2)? . Solving these 
equations yield nz, where n, = [(n, sin 64)? + aI. A similar derivation holds at each 


: ae . f , _211/2 . 

interface, resulting in the recurrence relationship n; = [Ci -sin ĝi-1)? + B; 2] í „i> 2.In 
addition, because the incidence angle of the laser is the same for each layer of the dielectric, 
6, =0, =: = 0, =m/2—a, the recurrence relationship can be simplified to n; = 


-311/2 -371/2 
[(nj_1 "sin 8;-1)? + 6°] e [(nj_1 + cosa)? + 6°] / ,i > 2. In summary, the refractive 
index and thickness of each layer of the dielectric, starting from the second layer, can be derived 
as a function of the parameters of the first layer. 


. = [(n;_1 + cos a)? + B; 
Li = L, X (tana)*/(cos a)i- 

In the first row of Equation (2), n; is the refractive index of the layer, and the label, i, 
follows the convention in Figs. 1 and 2. f; is the theoretical design velocity of electrons in the 
vacuum acceleration region on Surface I below the ith layer of the dielectric, which is also the 
velocity of light in the x-direction when it is in this dielectric. For this theoretical electron 
velocity, as the electron beam stream advances longitudinally in the acceleration channel, the 
electron energy and velocity increase. Electron velocity v, gradually increases when the 
electron moves below the surface of each dielectric layer, whereas the velocity v, of the 
corresponding laser pulse remains unchanged. Therefore, the phase difference gradually 
accumulates as the electrons move with the laser during the acceleration phase, ® = wt(1 — 
nf sin a). If the electron beam is injected at the maximum amplitude of the laser, it becomes 
necessary to ensure that the phase difference |®| < 7/2 when the electron beam leaves; 
otherwise, it enters the deceleration phase. This corresponds to two cases. (1) The electron 
velocity ve is extremely high, resulting in phase difference ® < — 7/2, and the electron 
beam reaches the front deceleration phase. (2) The electron velocity v, is extremely low, 
resulting in phase difference ® > 7/2, and the electron beam will enter the rear deceleration 
phase. 


-27/2 
i22 (2) 
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A suitable electron beam velocity, as expressed by Equation (3), can be selected to extend 
the acceleration distance and duration. This is achieved such that (1) the electron beam velocity 
is lower than the speed of light in the dielectric at the entrance of each dielectric layer, and (2) 
the electron beam velocity is higher than the speed of light in the dielectric at the exit of each 
dielectric layer. In Equation (3), x, is the longitudinal position of the electron, n is the entry 
into the nth layer of the dielectric, and YiijL;=0 when n=1. Thus, the relative 
longitudinal displacement of the electron beam stream in the phase is first backward 
(equivalently, the negative x-direction) and then forward (equivalently, the positive x-direction), 
extending the duration in the accelerated phase [9,15]. 


n-1 
Ve < Vo, Xe = } Li 
i=1 


I g 
Ve > Vo, Xe =) Hi 
i=1 


In the second row of Equation (2), L; is the thickness of the layers, known as Lı and L3 
in Fig. 1. From recursive formula (2), as long as the prism base angle a, refractive index n4, 
and thickness L, of the first layer of the dielectric are specified, the remaining refractive index 
n; and thickness L; can be obtained recursively. Moreover, Equation (1) indicates that when 
the incident electron energy is determined, only one degree of freedom remains between a and 
nı. Generally, researchers tend to select quartz and other common materials within the range 
of refractive indices and make the base angle æ as small as possible. A positive correlation 
exists between the acceleration electric field or acceleration gradient and cosa; hence, a 
smaller base angle a indicates a higher acceleration gradient [17-19]. In addition, the 
continuous or arbitrary regulation of the dielectric refractive index has led to significant 
progress in research on optical materials; therefore, a dielectric satisfying Equation (2) can be 
fabricated [21-23]. 


(3) 
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Fig. 2. Schematic of structure of 50 keV electron acceleration. The light blue, purple, and red ellipses represent the 
injected low-energy (50 keV) electron beam stream, accelerating electron beam stream (schematic), and output high- 
energy electron beam stream, respectively. The insertions “dec” and “acc,” indicate that half of the phases in a laser 
cycle accelerate the electrons, and the other half decelerate them. 
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Based on the theoretical analysis above, a two-dimensional numerical simulation study was 
conducted using the finite element method and the particle-in-cell algorithm. The simulated 
structure is shown in Fig. 2. The base angle a of the prism was set to 30°, and the refractive 
index and thickness of each layer were calculated using Equations (2) and (3), respectively. A 
symmetric double prism was used, in which the transverse forces on both sides of the channel 
centerline are of equal magnitude but in opposite directions. This configuration prevents beam 
deflections. The geometric parameters of the simulated structures are listed in Table 1. The 
laser was set as a-10.00 — um far-infrared CO2 laser. Because the timescale of the simulated 
physical process is similar to the pulse width of widely used CO2 lasers [24-28], the laser was 
set as a continuous plane wave with an electric field amplitude of Ey = 5.0 GV: m~t. The 
selected value implies that the power density of the laser beam is approximately 6.6 x 
101? W/cm?, and the energy density of the laser pulse irradiating the prism is approximately 
0.2 J/cm?. These values are consistent with the laser parameters used in a related study [17]. 
We propose that line-shaped laser pulses collimated by a cylindrical lens be utilized to drive 
the symmetric prism, enabling long-distance synchronization. With an energy of 50.00 keV, 
the injected electrons were in the low-energy nonrelative range (B « 1) typical for 
monoenergetic electron beams generated from electron guns with no initial energy dispersion 
and a charge of 9.36 fC. The particles in the bunch followed a KV distribution. The KV 
distribution uniformly places the particles in the phase space. The distance between the particles 
was approximately the same throughout the beam. The initially normalized emittance was 
103 pm: rad. The initial maximum transverse distance relative to the center of the beam was 
0.1 um, and the maximum divergence angle was approximately 23 mrad. The longitudinal 
beam length was 3.34 fs, and the distance between the center of the electron beam and prism 
Surface I was 0.5 um. These design parameters are consistent with those of the electron beam 
generated using the electron gun of the DLA [29-31], and the covered parametric range has not 


yet been investigated for the inverse Cherenkov-effect particle accelerator [16-19]. 
Table 1. Simulation parameters of geometric structure 


s refractive 
section length (um) ‘index 
first 31.31 3.96 
second 10.44 3.74 
third 13.92 3.56 
fourth 18.56 3.40 
fifth 24.74 3.21 
Structural 1.00 1.00 
spacing 


3. Simulation results and discussion 

The energy spectrum of the accelerated electron beam stream is shown in Fig. 3(a). Fig. 3(a) 
corresponds to the instant at 600.42 fs. Here, the electron beam had completely accelerated. 
The peak energy was 309.55 keV, and the full-width at half-maximum (FWHM) of the energy 
dispersion was 2.00 keV, corresponding to a relative ratio of 0.65%. Correspondingly, an 
electron energy gain of 259.55 keV was achieved with an average acceleration gradient of 
2.62 GeV: m7? andan acceleration factor Eacc/Eo ~ 0.52. This was obtained by dividing the 
total energy gain of the electrons by the entire path length of approximately 98.96 um. The 
phase-space diagram of the horizontal (y-direction) space is shown in Fig. 3(b), which 
corresponds to the same moment as in Fig. 3(a). The maximum divergence angle of the electron 
beam did not exceed 21 mrad. Compared with the initial value of approximately 23 mrad, 
the divergence angle was reduced, and a weak focusing effect was achieved. Fig. 3(c) shows 
the spatial (x-direction) of the electron beam and the energy distribution with time, with the 
sampling points separated by 0.257) or approximately 8.34 fs The electron energy 
increased monotonically, and continuous acceleration was achieved. The curve representing 
the electron beam energy increased smoothly through the vacuum channel, with only some 
jitter at the layer interfaces. This smooth increase indicates that the acceleration gradient does 
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not change significantly at each stage; thus, the proposed stacked acceleration structure 
achieves long-range acceleration with a high acceleration gradient. The electron beam 
underwent diffusion and recompression in the spatial x-direction, and the energy of the electron 
beam underwent diffusion and recompression. The stretched and diffused of the electron beam 
are caused by inhomogeneous field distribution in the accelerated vacuum channel, as reported 
in the literature [17-19]. When the diffusion of the electron beam stream becomes very strong, 
the positive phase (E, > 0), that is, the decelerating phase-electric field in front and behind the 
electron beam, suppresses it, introducing negative feedback. Consequently, the length of the 
accelerated electron beam is limited to half of the laser cycle, and the extremely-fast electrons 
slow down, inhibiting their energy increment and minimizing the energy dispersion. This also 
reduces the speed of the fast-electron beams and concentrates the electrons in the movement 
direction. This technique facilitates the electron beam to achieve its quasi-monoenergy and 
compression in the direction of motion after passing through the channel. This property may be 
valuable in future ultrashort pulse electron beam streams such as attosecond pulses, ultrafast 
electron diffraction, and other fields. In addition, the electron charge for the fluctuation of the 
electron beam current charge was observed (Fig. 3(d)), and the simulation demonstrated that 
no beam loss occurred during the electron beam current acceleration. This finding is consistent 
with the earlier analysis that the field generated by this stacked structure does not cause 
transverse (y-direction) deflection of the electron beam current. Furthermore, the width of such 
a vacuum channel is adequate. 
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Fig. 3. (a) Electron energy spectrum at end of acceleration, normalized to peak. (b) Phase-space diagram of 
horizontal (y-direction) space. The vertical coordinate is the divergence angle, and the horizontal coordinate is the 
lateral position of the electron. (c) Spatial (x-direction) and energy distribution of electron beam at each moment 
(difference of 0.25 Ty for each data sampling point). (d) Total power variation in electron beam at each moment 

(difference of 0.25 Ty for each data sampling point). 
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The spatial and temporal variations in the longitudinal electric field E, along the central 
axis of the vacuum channel were investigated to further examine the acceleration process of the 
electron. The path of an electron carrying average energy was adopted to approximate the mean 
trajectory of the electron beam, as shown in Fig. 4. This demonstrates how an electric field 
induces electron acceleration or deceleration. The blue-violet region represents the deceleration 
field, whereas the orange-red zone represents the acceleration field. The net energy gain of a 
particle is closely related to the phase in which the electron beam is placed in the electric field, 
E,. In the phase of the laser electric field, the solid green line represents the global line of the 
particle beam with average velocity. The electron beam is synchronized with the laser pulse 
and is always in the negative phase of the electric field (E,, < 0), known as the acceleration 
phase. This is why a stacked structure can achieve a high acceleration gradient for long-distance 
acceleration. Accelerating a single beam does not fully utilize all the acceleration phases. 
Therefore, in the future, a multibeam cluster acceleration mode similar to modern conventional 
accelerators, with each pair of neighboring electron beam clusters spaced by a deceleration 
phase, can be adopted to improve the energy conversion efficiency and electron beam power. 
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Fig. 4. Electric field E, (magnitude normalized with respect to Eg) as functions of time and space (x-position). The 
vector of this field is in the lateral direction of the vacuum channel. The green line indicates the electron trajectory, 
and the sampling time step is 0.257). The color bar on the right represents the normalized field strength, E,/Ep 
ratio. 


Moreover, the performance of each layer of the structure was assessed to adapt it to future 
chip accelerator applications and construct a structure that can be dismantled to alter the 
electron beam energy to satisfy the needs of various situations. The energy spectrum of the 
electron beam stream produced by the acceleration of each layer is shown in Fig. 5. The 
moments indicated by different colors in Fig. 5 correspond to the electron energy spectra 
obtained from the first layer to the sixth layer of the dielectric acceleration. The values are listed 
in Table 2. The energy dispersion of the intermediate electron beam flow was relatively high. 
However, none of the energy dispersions exceeded 4.00 keV of the FWHM, and none of the 
corresponding relative ratio values exceeded 2.00%. The peak energy of the electron beam 
obtained from the final acceleration was 309.55 keV, and the FWHM of the energy dispersion 
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did not exceed 2.00 keV. This indicates that the intermediate electron beam flow can also be 
used as the injection step, as discussed in the next section. 
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Fig. 5. Electron beam energy spectrum corresponding to end of acceleration of each layer of dielectric. 


Table 2 presents the results of an investigation of the acceleration efficiency, in which the 
energy gain in the electron beam stream in each accelerated stage is determined, together with 
the corresponding average acceleration gradient. This coincided with the case presented in Fig. 


3(c). 


Table 2. Electron beam profile and acceleration performance for each layer 


Average 
; electron half-height- i iy electron gain acceleration 
amon energy (keV) width (keV) ratio (%) energy (keV) gradient 
(GeV: m-t) 
first 125.55 0.34 0.27 75.55 2.41 
second 137.15 2.40 1.75 11.60 1.11 
third 178.75 3.53 1.97 41.60 2.99 
fourth 225.85 1.67 0.74 47.10 2.54 
fifth 309.55 2.00 0.65 83.70 3.38 
total 309.55 2.00 0.65 259.55 2.62 


This section presents the expandability of the proposed structure. First, the acceleration 
distance required to accelerate a 50 keV electron beam to 1MeV was approximately 
362.60 um using an average acceleration gradient of 2.62 GeV : m7? for the entire structure. 
Second, the thickness and refractive index of the dielectric to be added could be obtained using 
Equations (2) and (3). The length (x-direction or Surface I direction) and height (y-direction or 
Surface III direction, the highest point of the prism) of each subsequent layer are listed in Table 
3. In this case, the acceleration structure of the first five segments used the above structure. 
Furthermore, preliminary estimates were obtained for longer acceleration distances and higher 
energy gains. The electron relativistic velocity 6 and dielectric refractive index n changed 
with increasing electron energy (Fig. 6(a)). The pinch angle was based on the same initial angle 
as used in this study. The refractive index of the dielectric rapidly decreased as the electron 
energy increased, and the refractive index eventually reached 2.00 and remained constant. Fig. 
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6(b) shows the variation in the thickness of each section of the dielectric structure.The thickness 
of each stage increased exponentially as the number increased. Dielectrics of this size are easy 
to prepare. In addition, the laser could be infinitely extended in one direction (with the focal 
length or extendable distance significantly longer than the acceleration distance) and tightly 
focused in the other direction (laser wavelength magnitude) using beam spreading and line- 
shaped collimation [32]. This implies that on Surfaces I and II (Fig. 1), the passing distance of 
the laser in the dielectric is significantly longer than the required acceleration distance. In 
general, the proposed structure can be extended to higher energy targets by adding a new 


dielectric layer behind the structure. 
Table 3. Size of subsequent structure 


section length(um) height(um) 
sixth 32.99 76.18 
seventh 43.98 101.58 
eighth 58.64 135.44 
ninth 78.19 180.58 
tenth 104.26 240.78 
(a) == (b) ne 


Fig. 6. (a) Electron relativistic velocity (blue line) and dielectric refractive index (red line) as functions of electron 
energy; (b) Length of each section of dielectric structure (The green and brown lines indicate measurements in 
microns and meters, respectively, while logarithmic coordinates are used.). 


4. Conclusion 


A laser dielectric accelerator is proposed based on the inverse Cherenkov phenomenon using a 
multilayer structure. A parametric design approach was applied to this structure and is 
demonstrated with a design example for accelerating electron beam currents ranging from 
50.00 to 309.55 keV. This suggests that the stackable structure requires a high-energy-gain 
cascade acceleration technique featuring a high acceleration gradient and an extended 
acceleration distance, as demonstrated via the simulation. Furthermore, this stackable 
accelerating structure has a flexible configuration. Theoretically, the structure can achieve 
arbitrarily high acceleration under design assumptions by extending the accelerator length 
following the proposed design procedure. The proposed structure is a promising prototype for 
use in future chip accelerators. 
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